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A Non-potential Model for the Sun's Open Magnetic Flux 
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Abstract. Measurements of the interplanetary magnetic field (IMF) over several so- 
lar cycles do not agree with computed values of open magnetic flux from potential field 
extrapolations. The discrepancy becomes greater around solar maximum in each cycle, 
when the IMF can be twice as strong as predicted by the potential field model. Here we 
demonstrate that this discrepancy may be resolved by allowing for electric currents in 
the low corona (below 2.5Rq). We present a quasi-static numerical model of the large- 
scale coronal magnetic evolution, which systematically produces these currents through 
flux emergence and shearing by surface motions. The open flux is increased by 75% — 
85% at solar maximum, but only 25% at solar minimum, bringing it in line with esti- 
mates from IMF measurements. The additional open flux in the non-potential model arises 
through inflation of the magnetic field by electric currents, with super-imposed fluctu- 
ations due to coronal mass ejections. The latter are modelled by the self-consistent ejec- 
tion of twisted magnetic fiux ropes. 
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1. Introduction 

The "open magnetic flux" is the part of the Sun's mag- 
netic field that is attached at one end to the solar photo- 
sphere but extends out to form the interplanetary magnetic 
field (IMF). It is of practical interest because it drives vari- 
ations in the IMF near Earth. Short-term fiuctuations are 
a source of geomagnetic activity [Pulkkinen, 2007], while 
longer-term, secular variations — in particular, an apparent 
doubling over the last 100 years [Lockwood et al, 1999] — 
modulate the galactic cosmic ray fiux [Beer, 2000] and are 
recorded in terrestrial records such as tree rings or ice cores. 
It has even been suggested that these variations may be 
linked to changes in global climate through variation of cloud 
cover [Svensmark and Friis-Christensen, 1997]. 

Empirical models are able to reproduce the long term 
variation of the open fiux using either geomagnetic [Lock- 
wood et at, 1999, 2009b] or sunspot [Solanki et al, 2000; 
Vieira and Solanki, 2010] data as input. However, under- 
standing the origin of the Sun's open magnetic fiux requires 
knowledge of the structure and evolution of the magnetic 
field in the low solar corona. Unfortunately, direct magnetic 
measurements are available only at the photospheric level, so 
theoretical models are needed to extrapolate the magnetic 
field in the corona. The most popular approach has been 
to assume a potential field, i.e., a vanishing electric cur- 
rent density. This gives a unique magnetic field solution in 
the region between an observed radial magnetic distribution 
on the photosphere r = Rq and an outer "source surface" 
r ~ Rss where the field becomes purely radial. It is the 
magnetic flux through this source surface that determines 
the "open fiux" in the model. Based on comparison of PFSS 
models with observed coronal structures, it is conventional 
to set Rss = 2.5Rq [Schatten et al., 1969; Altschuler and 
Newkirk, 1969], although this leads to one criticism of the 
PFSS model: the real Alfven radius, above which the solar 
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wind kinetic energy density dominates the magnetic energy 
density, is of the order 10-Rq [Marsch and Richter, 1984]. 
Nevertheless, PFSS extrapolations with Rss = 2.5i?Q are 
able to reproduce the locations of many He 10830A coronal 
hole boundaries [Wang et al, 1996; Schrijver and DeRosa, 
2003], which map the distribution of at least some open field 
line footpoints on the solar surface. However, when the mag- 
nitude of open fiux in such PFSS extrapolations is compared 
with in situ IMF measurements, it tends to be too low. 

To compare the predicted open fiux with in situ observa- 
tions of the IMF strength, the modelled open fiux is mapped 
out to 1 AU, and it is assumed that non-radial motions in 
the intervening space act to even out the distribution in lat- 
itude so that it becomes uniform [e.g., Wang et al, 2000]. 
This latter assumption is based on out-of-ecliptic magnetic 
field measurements by Ulysses, which show that the radial 
heliospheric flux is essentially independent of heliographic 
latitude [Balogh et al., 1995; Smith et al., 2001]. The pre- 
dicted radial field strength from the model at 1 AU, , is 
then given simply by scaling the total unsigned flux through 
the outer model boundary, i.e. 



A-kR% 



[Br{Rss,e,<t>)\dn, 



(1) 



where Re the radius of Earth's orbit. Several authors have 
used observed photospheric magnetograms to produce time 
series of PFSS extrapolations, comparing B^ with the ob- 
served Bx component of the IMF [Wang and Sheeley, 1995; 
Wang et al, 2000; Schrijver and DeRosa, 2003; Schiissler 
and Baumann, 2006]. In Figure 1, the grey line shows the 
observed Bx for Solar Cycle 23, while the three colored lines 
show PFSS predictions of B^ using Equation 1 with differ- 
ent photospheric magnetic fleld data (Section 2). Section 2 
discusses these PFSS models in more detail; the important 
point is that all three PFSS predictions are too low com- 
pared to the observations, particularly during more active 
phases of the Cycle. 

This shortfall has been noted over the past three solar 
cycles [e.g., Riley, 2007], and several possible explanations 
have been proposed. Firstly, the open flux may be increased 
simply by lowering the source surface -Rss in the PFSS 
model below 2.5Rq, but this leads to poor agreement with 
coronal structures observed during eclipses [Altschuler and 
Newkirk, 1969], and Rss should arguably be larger than 
2.5-Rq, not smaller [Jiang et al, 2010]. It is possible that 
the source surface is not spherical and its shape varies over 
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Figure 1. Comparison of measured radial IMF at 
lAU with various PFSS extrapolations over Solar Cy- 
cle 23. The thick grey line shows daily 0MNI2 \B^\ 
data with a two-stage averaging [see Lockwood et ai, 
2009b]: an initial daily average of the signed data, to 
smooth out local small-scale fluctuations, followed by a 
27-day running Carrington average of the unsigned data. 
Coloured lines show PFSS predictions taken from WSO 
with latitude-dependent correction (purple), WSO with 
constant correction (red), and NSO/Kitt Peak (green). 
The smoothed monthly sunspot number (from SIDC) is 
shown in blue, with letters A to D indicating times of our 
NP simulations. 

the cycle [Schulz et al., 1978; Riley, 2006], but the effect of 
this on the open flux is presently unclear. Secondly, a modi- 
fied version of the PFSS model, the so-called "current-sheet 
source surface" (CSSS) model, is found to give better agree- 
ment with the observed IMF [Zhao and Hoeksema, 1995; 
Schiissler and Baumann, 2006]. The CSSS model adds an 
inner spherical "cusp surface" r = Res beneath the outer 
source surface (which is moved out to about IORq). Beneath 
Res (typically about 1.8Rq) the field remains potential, 
but currents are introduced between Res and Rss- This 
technique reproduces the latitude-independent field found 
by Ulysses as well as the heliospheric current sheet. How- 
ever, the CSSS model has essentially three free parameters 
(Res, Rss, and a parameter controlling the current distri- 
bution), which may be varied in order to match the open 
flux to the observed IMF [e.g., Jiang et ai, 2010]. A third 
possible explanation for the PFSS shortfall in open flux, pro- 
posed by Fisk and Zurbuchen [2006], is that reconnection of 
open field lines with closed coronal loops could lead to a 
"diffusion" of open flux into otherwise closed regions, pro- 
ducing an additional component of open flux not included 
in PFSS models. However, this has not been demonstrated 
by observations [Rust et al., 2008]. A fourth possibility, re- 
cently suggested by Lockwood et al. [2009a, b], is that flow 
speed variations in the solar wind above 2.5Rq could lead 
to an increase in the unsigned radial flux at increasing dis- 
tance from the Sun. Such a "flux excess" effect appears to 
be supported by observations from various spacecraft out 
to 20 AU, and could potentially also explain the difference 
between PFSS flux at 2.5Rq and IMF measurements. 

In this paper, we suggest another possible explanation for 
the PFSS shortfall: the influence of electric currents in the 
low corona below 2.5Rq. Such currents are expressly pro- 
hibited in the PFSS model, yet we may infer their presence 



throughout the solar cycle, both in newly-emerged active 
regions, e.g., in X-ray sigmoids [Canfield et ai, 1999, 2007], 
and outside active latitudes, including long-lived Hq fila- 
ment channels [Martin et al, 1994; Mackay et ai, 2008] and 
coronal magnetic flux ropes [Gibson and Fan, 2006]. There 
have already been indications from global fuU-MHD mod- 
els that allowing such currents will increase the open flux 
[Riley, 2006]. However, these models tend to assume an 
initial potential field which is then relaxed to equilibrium 
with an imposed solar wind structure, so do not take into 
account the time-dependent development of currents due to 
the emergence and interaction of active regions in the lower 
corona. In Section 3, we describe an alternative, simpli- 
fied global magnetic field model that is less computationally 
demanding yet takes these time dependent non-potential ef- 
fects into account. The remainder of the paper applies this 
model to study the magnitude and variation of the open 
flux. In Section 4 we consider the effect of the outer bound- 
ary condition on the open flux, before showing in Section 5 
how the NP model predicts enhancement of the open flux 
due to inflation of the magnetic field by coronal currents. 
Conclusions are given in Section 6. 

2. Open Flux in the PFSS Model 

With a fixed source surface radius Rss, the open flux in a 
PFSS extrapolation depends purely on the photospheric ra- 
dial magnetic field at a single instant. This is the cumulative 
product of many bipolar active regions which have emerged 
and been dispersed by surface flux transport. In fact, the 
lowest order multipole components of the photospheric field 
determine the open flux, since higher-order multipoles decay 
more rapidly with height [Wang and Sheeley, 2002; Mackay 
et al, 2002]. For a few Carrington rotations, newly-emerged 
bipoles lead to fluctuations of the non-axisymmetric, equa- 
torial dipole component. Later, however, this contribution 
decays, or cancels between neighbouring bipoles, and the 
net asymptotic contribution is to the axisymmetric dipole 
component. During a single Solar Cycle, there is a factor 
2 modulation of the observed IMF {e.g., Figure 1), though 
the variation is typically only 10%-20% that of the closed 
flux [Wang et al., 2000]. The observed open flux is not quite 
in phase with sunspot activity (Figure 1), peaking about 1- 
2 years after Maximum. Earlier PFSS extrapolations from 
surface flux transport models had trouble reproducing this 
observed phase relation [Mackay and Lockwood, 2002], al- 
though this might be explained either by the properties of 
active regions used [Schiissler and Baumann, 2006]), or by 
including the contribution of CMEs, discussed in Section 
5.2. 

Because it is determined by the photospheric magne- 
togram input, the open flux predicted by PFSS models de- 
pends on the source of these observational data, leading to 
differing predictions using magnetograms from different ob- 
servatories. For observations made in the Fe 5250 A line at 
either Mount Wilson Observatory (MWO) or Wilcox Solar 
Observatory (WSO), a major factor affecting the level of 
open flux predicted is the saturation correction applied to 
the magnetograph signal. To illustrate this point, we have 
computed PFSS extrapolations from WSO data with two 
alternative correction factors that have been used in the lit- 
erature. The red line in Figure 1 uses a constant factor 
1.85, which is appropriate for the WSO instrument [see Ri- 
ley, 2007], whereas the purple line uses a latitude-dependent 
factor 4.5 — 2.5 sin^ A that was derived from MWO measure- 
ments of the same line [Ulrich, 1992; Wang and Sheeley, 
1995]. Clearly the latitude-dependent correction factor gives 
closer agreement to the observed B^ , even though it is ar- 
guably not applicable to the WSO instrument [Riley, 2007]. 
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For comparison, the green hne shows the PFSS prediction 
using alternative magnetogram data from NSO/Kitt Peak 
derived from the Fe 8688A line, where a similar saturation 
factor is not required [Arge et al., 2002]. This higher reso- 
lution data is used to derive the emerging bipole data input 
to our non-potential simulations. 

It is clear from Figure 1, however, that no matter which 
data are used, PFSS extrapolations predict too low an open 
flux over much of the cycle, particularly in more active 
epochs (such as the years 2000 or 2004). Agreement is better 
in Minimum periods, although more so in 1996 than in 2008. 
The non-potential enhancement described in this paper is 
one possible explanation for this solar activity-dependent 
discrepancy. 

3. Non-potential Evolution Model 

With the aim of moving beyond the PFSS model and 
understanding the evolution of currents in the coronal mag- 
netic field, we have developed global, non-potential simu- 
lations based on the coupled flux transport and magneto- 
frictional model of van Ballegooijen et al. [2000]. This non- 
potential (hereafter NP) model follows the time evolution 
of the large-scale magnetic field through a quasi-static se- 
quence of near force-free equilibria, in response to flux emer- 
gence and shearing by surface motions. Currents are gener- 
ated in the corona, and subsequent flux cancellation above 
polarity inversion lines tends to concentrate the associated 
magnetic helicity into either sheared arcades or twisted mag- 
netic flux ropes [through the mechanism described by van 
Ballegooijen and Martens, 1989]. The model in its present 
form has been successfully applied to study the formation 
of filament channels [Mackay et al, 2000; Yeates et al., 
2008a] and the initiation of coronal mass ejections [Yeates 
and Mackay, 2009]. The latter occur in the model when 
flux ropes grow too strong to remain in equilibrium [Forbes, 
2000] and are ejected through the outer boundary of the 
numerical domain. 

The NP model differs firom PFSS extrapolations in a 
number of important respects. Firstly, the latter assume 
a current-free corona, whereas in the NP model currents are 
generated both by the emergence of active regions and by 
subsequent shearing by photospheric motions. The effect 
of these currents on the open flux is to increase it relative 
to a potential field. Secondly, the PFSS model uses a se- 
quence of independent extrapolations, usually only once per 
solar rotation (27.27 days), with no information on how the 
coronal or photospheric fields have evolved. In contrast, 
the NP model follows the continuous evolution so that the 
magnetic field can retain a "memory" of previous interac- 
tions and topology. Because an observed magnetogram of 
the whole solar surface can be taken only once per rotation, 
we do not use synoptic magnetograms directly as input into 
the NP model, but rather use them to determine where on 
the solar surface new active regions have recently emerged. 
Corresponding magnetic bipoles axe then inserted individu- 
ally into the simulation to produce a continuously evolving 
boundary condition. A further difference is at the outer 
boundary. To facilitate direct comparison with PFSS ex- 
trapolations, this boundary is kept at r — 2.5Rq in this pa- 
per. However, the "source surface" condition used in PFSS 
extrapolations, namely that the transverse magnetic field 
components vanish, can not be apjjlied in the NP model 
since it would prevent the ejection of magnetic flux ropes 
through this boundary following their loss of equilibrium. 
Instead, we impose a radial outflow to simulate the effect 
of the solar wind in radially opening field lines that reach 
near to the boundary [Mackay and van Ballegooijen, 2006a]. 
This is discussed in Section 4. 



3.1. Model Equations 

The simulations use a domain extending from 0° to 360° 
in longitude, —80° to 80° in latitude, and Rq to 2.57?© in 
radius [Yeates et al, 2008a]. Note that, to avoid the compu- 
tational difficulties of grid convergence and coordinate sin- 
gularity, our domain omits the regions within 10° of the 
poles. These omitted regions amount to only 1.5% of the 
solar surface. Although they may carry a larger fraction 
of the open flux at solar minimum, this is partially offset 
by displacement of flux to lower latitudes, and may well be 
within the observational uncertainty of the magnetograms 
at these high latitudes. 

The large-scale mean magnetic field is evolved by the in- 
duction equation 



dAo 
dt 



: vo X Bo -I- £o, 



(2) 



where Ao(r,t) is the vector potential for the mean mag- 
netic field B() = V X A(), with gauge chosen to cancel the 
additional gradient term on the right-hand side of (2), and 
Vo (r, t) is the mean plasma velocity. The mean electromo- 
tive force £o{r,t) describes the effect of small-scale fiuctua- 
tions which are not resolved in our mean field model, such 
as braiding and current sheets produced by interaction with 
convective fiows in the photosphere. We assume the form 



So 



-mo, 



where jo = V x Bo is the current density, and 



?? = % 1 + 0.2 



iJol 
|Bo| 



(3) 



(4) 



is a turbulent diffusivity. The first term is a uniform back- 
ground value rjo = 45km'^s^^ and the second term is an 
enhancement in regions of strong current density [Mackay 
and van Ballegooijen, 2006a], introduced to limit the twist 
in helical fiux ropes to about one turn, as observed in fila- 
ments [Su et al, 2009]. 

Rather than solving the full MHD equations for the ve- 
locity vo, we approximate the momentum equation in the 
coronal volume by the magneto- frictional method [Yang et 
al., 1986], setting 



1 jo X Bo 
^° = M — • 



(5) 



This artificial velocity models the coronal field evolution 
through a sequence of near (nonlinear) force-free equilibria. 
In addition, a term is added to the radial component of vq 
near the outer boundary; this will be described in Section 4. 

On the lower, photospheric boundary, the magneto- 
frictional velocity is not applied, but rather the radial mag- 
netic field, Bor, is evolved using the standard surface fiux 
transport model [Sheeley, 2005, and references therein]. 
This surface component of our simulation was described in 
detail in Yeates et al. [2007]. In terms of the vector poten- 
tial Ao and standard spherical polar coordinates {r,6,4>), 
the evolution equations on the lower boundary r = Rq are 



dAn 



dt 



= U^Bor 



= —UgBor + 



D dBor 
Rq sin 9 d(j) 
D dBor 



Rq 



■0 



(6) 
(7) 



Here D = 450km'^s~^ is a diffusivity modelling the ran- 
dom walk of magnetic fiux owing to the changing super- 
granular convection pattern. The differential rotation veloc- 
ity = Q,{6)RQS\n6 uses the observationally-determined 
Snodgrass [1983] profile 



fi(^) = 0.18 - 2.3 cos^ e - 1.62COS* 6ldegday"\ 



(8) 
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Figure 2. Comparison of tiie NP model and PFSS extrapolations, for periods A to D. The plots on 
the left show the predicted and observed \Bx\, while the snapshots show the magnetic field in the 
NP model (middle) alongside a PFSS extrapolation (right) from the same radial photospheric field. In 
the left-hand plots, as in Figure 1, grey curves show the observed \Bx\, thick black curves show pre- 
dictions from the NP simulations, the black dashed line shows the PFSS prediction from the simulated 
photospheric field, and coloured symbols joined with dashed lines show PFSS predictions from observed 
magnetograms. The vertical dotted lines indicate the approximate end of the initial "ramp-up" period. 
In the snapshots, grey-shading shows radial magnetic field on the photosphere (white positive, black 
negative, with a saturation level of 20 G or 2mT), while selected coronal magnetic field lines are traced 
from the plane of the sky. 



written in the Carrington frame. We include a poleward 



as in Mackay and van Ballegootjen [2006a]. The constant C 
is chosen to give a peak flow at mid-latitudes of 16ms~^, 
and the flow vanishes at the boundaries ^jjjjjj = 10° and 
Smax = 170° of the computational domain. 



bulk flow, or meridional circulation 



ue{6) — C cos 



max + 



261) 



2 (emax 6ijjjin) 



cos 6, 



(9) 



3.2. Observational Input 

Newly emerging magnetic bipoles are inserted based on 
active regions observed in synoptic normal-component mag- 
netograms from the US National Solar Observatory (NSO) 
at Kitt Peak. The locations of these regions on the solar sur- 
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Table 1. Simulation periods. 



Period 


Carrington Rotations 


Dates 


Bipolcs 


Instrument 


A (Minimum) 


1911.5 to 1917 


12 Jul 1996 to 04 Jan 1997 


16 


KPVT 


B (Maximum) 


1962.5 to 1968 


03 May 2000 to 27 Oct 2000 


122 


KPVT 


C (Declining phase) 


2018.5 to 2024 


08 Jul 2004 to 02 Jan 2005 


52 


SOUS 


D (Minimum) 


2067.5 to 2073 


06 Mar 2008 to 30 Aug 2008 


14 


SOLIS 



face are determined from the magnetograms, but the exact 
date of emergence cannot be determined from the available 
observations, since at any time these show only the earth- 
ward face of the Sun. We simply insert all new regions 7 
days before their date of first observation at central merid- 
ian, with their properties appropriately scaled [as in Yeates 
et al., 2007]. Varying this time period does not change the 
accuracy of the surface magnetic field, and does not affect 
the level of the open flux on timcscales of a solar rotation 
or longer [Schrijver and DeRosa, 2003]. The insertion is 
instantaneous, but is preceded by an artificial "sweeping" 
of strong magnetic field out of the insertion region, as de- 
scribed in Yeates et al. [2008a]. This models the expected 
distortion of preexisting coronal field by a newly-emerging 
flux tube, and prevents the formation of disconnected flux 
in the corona. We find that this technique of inserting indi- 
vidual magnetic bipoles to model emerging active regions is 
able to retain high accuracy between the sirrmlated and ob- 
served surface fields over many months [Yeates et al., 2007]. 

The emerging bipolcs take the mathematical form given 
in Yeates et al. [2008a, equations 6 to 9], with properties 
chosen to match the size, tilt, and magnetic flux of the ob- 
served regions. They arc inserted in 3D, and are given a 
non-zero twist (magnetic helicity) through the parameter /?, 
as described in Yeates et al. [2008a]. The magnitude and, 
especially, sign of this emerging bipole helicity have been 
shown to influence both the chirality of filament channels 
[Mackay and van Ballegooijen, 2005; Yeates et al., 2008a] 
and the ejection rate of magnetic fiux ropes [Yeates and 
Mackay, 2009]. Unfortunately, the optimum value of the 
helicity parameter for each bipole is poorly constrained by 
present-day observations, and must be arbitrarily chosen. 
The observations that do exist indicate a variation in both 
magnitude and sign of helicity both within individual active 
regions and between different regions, as well as a nega^ 
five gradient of helicity with latitude [Pevtsov et al., 1995]. 
To approximate these features, we select a random /? value 
for each inserted bipole (up to a maximum j^Sj = 1) from 
a normal distribution with mean — 0.4A2/25° and standard 
deviation 0.4, where Ao is the central latitude of the bipole. 

3.3. Four Simulations over Cycle 23 

In this paper we focus on four main simulation runs A 
to D, each covering 5.5 months at different epochs of Cy- 
cle 23. The dates are shown in Table 1, and the periods 
are indicated in Figure 1. The fourth column of Table 1 
shows the total number of bipolcs emerged in each simula- 
tion run, as determined from the observed magnetograms. 
The right-most column indicates the instrument used; the 
older vacuum telescope (KPVT) was replaced in 2003 by SO- 
LIS (Synoptic Optical Long-term Investigations of the Sun). 
The four runs A to D arc illustrated in Figure 2, where plots 
in the left column compare the modelled from both the 
PFSS and NP models with the observed IMF jB^^^j. The 
snapshots in the middle and right columns show the NP 
simulation (left) and a PFSS extrapolation from the same, 
simulated, surface radial field (right). The black dashed lines 
in the left column of Figure 2 show for daily sequences of 
such PFSS extrapolations. Their agreement with the PFSS 
extrapolations from observed magnetograms (once per solar 
rotation; colored dashed lines) indicates that our technique 
for incorporating observed bipoles and simulating the sur- 
face field is sufficiently accurate to study the open flux. 

It is evident firom Figure 2 that the NP model for each 
period, after an initial ramp-up over about 1.5 months (in- 



dicated by the vertical dotted lines) , settles on a level of B^ 
which more closely matches IMF observations than that ob- 
tained from PFSS extrapolations. The initial ramp-up (to 
the left of the vertical dashed lines) demonstrates the im- 
portance of time evolution. This is the time taken for newly 
emerging bipoles and surface motions to evolve the coronal 
field away from the initial condition (a PFSS extrapolation) 
to some sort of equilibrium level of non-potentiality. To un- 
derstand how the NP model achieves the observed level of 
open fiux we must consider several contributory effects. 

4. Effect of the Outer Boundary Condition 

In our NP model, a radial outflow is imposed at the outer 
boundary r = 2.5Rq to model the effect of the solar wind 
in opening up coronal field lines above this height. A term 
of the form 

vout = ^'oc-(^««o-'-)/'-"f (10) 

is added to Equation 5, where the exponential fall ofi' (width 
r^i = 0.2Rq) ensures that only structures near to the outer 
boundary are affected. Typically wo set vo = lOOkms"^ 
as a compromise between a realistic solar wind speed and a 
reasonable time step in the integration (at least 30 s). 

Because it affects only those closed fields already near to 
the outer boundary, the speed vq has a relatively small effect 
on the open fiux in the NP model. Physically, this corre- 
sponds to the plasma-/3 at lower heights being too low for 




1999.3 1999.4 1999.5 1999.6 1999.7 
Year 



Figure 3. Effect of a radial outflow at the outer 
model boundary, illustrated for a period in 1999. Thick 
lines show the predicted B^ in simulations with vo = 
lOOkms"^, while thin lines show that for vo = 50kms~'^. 
The two solid black lines show the NP simulations, the 
two dot-dashed lines show PFSS extrapolations relaxed 
to equilibrium with the imposed outflow, and the dashed 
line shows the PFSS prediction in the absence of any out- 
flow. Note that all PFSS extrapolations in this plot are 
taken from the simulated photospheric field. The thick 
grey line shows the observed daily OMNI2 \Bx \ data with 
a 27-day running average. 
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the plasma flow to deform the magnetic field. This small ef- 
fect is demonstrated in Figure 3, which shows predicted 
from various coronal models using the same photospheric 
field for 5 months in 1999 (along with the observed Bx in 
grey). The two solid black lines show runs of the NP sim- 
ulation with different outfiow speeds, vo ~ 50kms~^ (thin 
line) and lOOkms"^ (thick line). Halving the outfiow ve- 
locity is seen only to decrease the level of open fiux by 
about 10%. The difference is visible in equilibrium levels 
are reached in the first few days of the simulation, before 
any bipoles emerge (the initial conditions are a PFSS ex- 
trapolation, hence the initial increase). Another effect seen 
in Figure 3 is a right-ward shift of the jagged peaks in open 
fiux, which correspond to flux rope ejections, in the run with 
lower outflow speed. This is because the outfiow affects the 
evolution of fiux ropes after they lose equilibrium (though 
not before). With lower outfiow they take longer to pass 
through the outer boundary. 

To further isolate the effect of the outflow, we have taken 
PFSS extrapolations at several times, inserted these as ini- 
tial conditions in the non-potential code, and relaxed to 
equilibrium with an imposed outflow (but no flux emer- 
gence or surface transport). Two such sequences, with 
Vo = 50kms~^ and wo ~ lOOkms"^, are shown by the 
dot-dashed lines in Figure 3. While there is an enhance- 
ment of 30%-40% relative to the original PFSS extrapola- 
tion (dashed line), the open flux is still much lower than in 
the full NP simulation, which better matches the observed 
IMF. The radial outflow is therefore not responsible for most 
of this improved agreement. Other effects are important and 
are considered in the next section. 

5. Open Flux in the Non-potential Model 

The total open flux in the NP model originates from three 
main sources, in addition to a small dependence on the out- 
flow boundary condition as described in Section 4: 

1. The photospheric flux distribution, which originates in 
the observational magnetogram input. 



2. A long-term enhancement to the mean level of open 
flux, due to inflation by coronal currents. 

3. Fluctuations super-imposed on the mean enhance- 
ment, caused by the ejection of magnetic flux ropes. 

The first source also applies to PFSS extrapolations and 
was described in Section 2. Here we consider the other two 
sources in turn. 

5.1. Long-term Enhancement by Coronal Currents 

This is a consequence of the basic tendency of mag- 
netic flux tubes to expand as they are stressed by cur- 
rents [Mackay and Lockwood, 2002; Sturrock et at, 1994]. 
The early numerical force-free model of Barnes and Stur- 
rock [1972] showed that held lines expand outwards as the 
magnetic held is twisted up. In a comparison between 
global MHD and PFSS solutions, Rtley [2006] found that the 
magnetic held in the MHD solution was more inflated and 
opened, with a realistic "cusp-like" morphology of coronal 
streamers not present in the PFSS solution. Such features 
are also seen in our NP simulations (Figure 2). 

In our NP model there are three main sources of currents, 
in addition to the radial outflow discussed above. These 
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Figure 5. Photospheric footpoints of open field lines 
on 27 October 2000, in the NP model (blue) and PFSS 
extrapolation (from the simulated surface field, red out- 
lines). Grey shading shows the radial magnetic field in 
the photosphere (white positive, black negative). 
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Figure 4. The four main sources of coronal currents in 
our NP model. 
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Figure 6. Dependence of the ratio of open fiux on the 
ratio of NP to PFSS energy, for simulation periods A to 
D. Data points are plotted for every second day of each 
simulation. The ratio of open flux is defined as the ratio 
of total unsigned radial fiux through the source surface 
r — 2.5-Rq to that through the photospheric boundary 



YEATES ET AL.: THE SUN'S OPEN MAGNETIC FLUX 



X- 7 




Figure 7. Sequence showing the loss of equihbrium of a magnetic flux rope in the NP simulation (period 
B). The interval between panels (a) and (b) is 4 days, with subsequent panels at 2-day intervals. Grey- 
shading shows radial magnetic field on the photosphere (white positive, black negative, with a saturation 
level of 20 G or 2 mT) , while selected coronal magnetic field lines are traced from the plane of the sky. 



are shown schematically in Figure 4 [see also Yeates et ai, 
2008b] . Firstly, the new bipoles emerge with a non-zero twist 
(a). This is initially concentrated low down in the center of 
the bipole, shearing field lines that cross the central polarity 
inversion line. Secondly, when bipoles emerge they displace 
pre-existing fields and produce currents at the interfaces be- 
tween old and new flux systems (b). This can lead to the 
formation of observed "intermediate" filaments at these lo- 
cations [Gaizauskas et al., 1997; Wang and Muglach, 2007; 
Mackay et al, 2008]. Thirdly, over days to weeks, surface 
motions, primarily differential rotation, shear the coronal 
field, generating further currents (c). Supergranular diffu- 
sion converges magnetic flux and liclicity toward polarity in- 
version lines, leading to flux cancellation and the formation 
of twisted, current-carrying magnetic flux ropes [ Yeates and 
Mackay, 2009]. Figure 5 shows that these different sources 
of current in the NP model lead to additional open field 
lines at various locations on the solar surface. The lower 
latitude open regions that were present in the PFSS model 
(red outlines) are much extended and their shapes modified, 
while additional open regions are found in a number of ac- 
tive regions. However, shearing by differential rotation at 
higher latitudes has also led to extension of the polar coro- 
nal holes. We do not discuss these spatial differences in the 
distribution of open field lines any further in this paper, con- 
centrating instead on the global, integrated open flux over 
time. In particular, spatial comparison of the NP model 
with observed coronal hole boundaries is left to future work. 

Figure 6 demonstrates how the fraction of open flux (i.e., 
the ratio of unsigned open flux to total unsigned photo- 
spheric flux) depends on the ratio of total magnetic energy 
in the NP and PFSS models. This energy ratio is a proxy 
for the total current and measures the "non-potentiality" of 
the global field. Each symbol represents a single snapshot, 
and the shapes indicate the four simulation periods as in 
the legend. Note that each run was initialized with a PFSS 
extrapolation, so initially has am energy ratio of unity. In 
each simulation run, there is a clear trend for a higher ra- 
tio of open flux at times with greater free magnetic energy, 
consistent with the idea of inflation by the presence of cur- 
rents. Another clear trend in Figure 6 is the solar cycle 
dependence, which manifests itself in two ways: energy ra- 
tios are typically lower at Minimum, but the ratio of open 
flux to photospheric flux is higher. Since the free energy and 
current arise from the emergence and interaction of strong 
active regions, it is no surprise that there is less free energy 
at Minimum. Indeed, the total parallel current in the simu- 
lations increases from period A to period B by about a factor 
5.5. The higher fraction of open flux at Minimum is a well 
known result, arising because at Minimum the open flux de- 
creases much less than the total photospheric flux [ Wang et 
al., 2000], owing to the presence of the polar coronal holes. 
We stress, however, that the absolute amount of open flux 



(i.e., without normalizing by the photospheric flux) is higher 
at Maximum, in accordance with IMF observations, as is ev- 
ident in Figure 2. It is this activity-dependent enhancement 
of the magnitude of open flux that is the key difference be- 
tween our NP model and the PFSS model, responsible for 
the improved match with observations. 

5.2. Fluctuating Enhancement by Flux Rope Ejections 

In addition to the long-term inflation of the magnetic 
fleld by coronal currents, there are additional temporary 
enhancements of the NP open flux lasting for a few days. 
Many of these are caused by the ejection of flux ropes, 
which form by flux cancellation at polarity inversion lines 
and lose equilibrium if they gain too much axial flux rela- 
tive to the overlying field. An example is shown in Figure 7. 
Using the same magncto-frictional model as in the present 
paper, Mackay and van Ballegooijen [2006b] looked in detail 
at the flux rope that formed and lost equilibrium between 
two bipoles in a simplified configuration. They found that 
the open flux of the two bipoles increased temporarily by 
a factor ~ 2 during the ejection as closed field lines were 
opened and reconnected. 

An analysis of the fluctuations in one of our global sim- 
ulations (period B) is presented in Figure 8. Figure 8(a) 
shows the total B (thin line) with a 14^day running mean 
overlaid in the thick line. This mean component arises from 
the various longer-term contributions as discussed earlier. 
Figure 8(b) shows the (absolute) time rate of change of 
, where several sharp peaks are clearly visible. To re- 
late these fluctuations to flux rope ejections in the model, 
Figure 8(c) shows the number of "erupting points" on each 
day in the simulation. These are points on a computational 
grid which (a) are identified as part of a flux rope structure, 
and (b) have a radial velocity in the magneto-frictional code 
in excess of 0.5kms~^. Details of our automated techniques 
to detect flux ropes and their ejections in the global sim- 
ulations are given by Yeates and Mackay [2009]. We see 
from Figure 8(c) that the number of grid points involved in 
flux rope ejections has a sharply-peaked structure similar to 
that of \dB^/dt\. These peaks correspond to the times of 
flux rope ejections, though some ejections take place grad- 
ually over a number of days, and there may be several in 
progress at one time. There is a reasonable correspondence 
between peaks in Figure 8(b) and 8(c) for example around 
2000.54 — although this is by no means one-to-one, demon- 
strating the complex nature of the global model. Computing 
the cross-correlation between the number of erupting points 
and the fluctuating component of B^ (deflned as the total 
minus the mean shown in Figure 8a), we flnd a peak correlar 
tion of 0.35 at a time lag of 1 day, with no strong correlation 
at later times (Figure 8d) . This supports our assertion that 
flux ropes enhance the open flux while they are ejected. 
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Comparing the diflFcrcnt NP runs, wc find that ejections 
occur at the rate of 0.11, 0.87, 0.49 and 0.19 per day in the 
final 100 days of runs A, B, C, D respectively. We use only 
the last 100 days here to avoid the initial ranip-up in the 
ejection rate (and also the open fiux). This solar cycle vari- 
ation in the number of flux rope ejections leads to a similar 
variation in contributions to the open flux. This may be 
seen in Figure 2, where the predicted for periods B and 
C shows more short-term fluctuations than that for periods 
A and D. This is in part responsible for the enhancement of 
open flux at Maximum which brings the NP model closer to 
the observed than the PFSS model, although the back- 
ground inflation by currents is also important. Mackay and 
van Ballegooijen [2006b] estimated that if flux rope ejec- 
tions occur once every 20 days from each active region, and 
enhance the open flux for 2 days, then this fluctuating con- 
tribution can account for only about 10% of the open flux. 

Observationally, flux rope ejections are likely to be man- 
ifested as observed coronal mass ejections (CMEs), and in- 
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Figure 8. Relation between the open flux and flux rope 
lift-offs in the NP model, for simulation period B. Panel 
(a) shows the total B^ (thin line), with the "mean" value 
after smoothing with a 14^day running average superim- 
posed (thick line). Panel (b) shows the absolute time 
rate of change of the (total) B^ , while panel (c) shows 
the number of erupting flux rope points present on each 
day (see text). Panel (d) shows the cross-correlation be- 
tween (c) and the fluctuating component in (a) . 



deed CMEs have been linked to freshly opening regions of 
coronal magnetic flux [Luhmann et al., 1998]. Owens et al. 
[2008] recognized the temporary contribution of CMEs to 
the IMF, and proposed that the IMF is the combination 
of a constant open flux "floor" and transient contributions 
from Earth-directed CMEs. They found a strong correlation 
between the IMF strength and Carrington rotation averages 
of the observed CME rate, and went on to estimate that if 
each CME contributes ~ 10^^ Wb of magnetic flux to the 
IMF, it must do so for a period of 30 50 days in order to 
explain the observed cycle variation of the IMF. While this 
flux per ejection is comparable in order of magnitude to 
ejections in our NP model, the timescale for contribution is 
much longer. By adding the estimated CME contribution to 
the open flux predicted from the PFSS model, Riley [2007] 
was able to match the observed IMF with a more realistic 
timescale of about 2 days. These are only rough estimates, 
however. Our NP simulations represent the first step toward 
a time-dependent coronal model that sclf-consistcntly initi- 
ates CMEs and the associated open fiux enhancement, al- 
though predicting the enhancement by individual observed 
events is not yet possible. Primarily due to the low time 
resolution of the magnetogram input to the simulations, 
the rate of ejections in the model is only about 20% that 
of observed CMEs, and the model is unable to reproduce 
the rapid onset of multiple CMEs in individual active re- 
gions [Yeates et al., 2010]. One might initially expect that a 
more finely-detailed model that produced the observed rate 
of ejections would yield too high an open fiux. However, this 
may be balanced by two factors. Firstly, many of the ad- 
ditional CMEs recur in the same active region within a few 
hours to days, so might not generate significant additional 
open flux. Secondly, the timescale for open flux enhance- 
ment by each ejection, which is often a number of days in 
our simulations, could be too long. More frequent but more 
rapid ejections may not lead to additional open fiux beyond 
our model. Future work will need to consider the dynam- 
ical effects controlling the ejection of fiux ropes after they 
lose equilibrium, which are not taken into account in the 
magneto-frictional model. 

6. Conclusion 

We have studied the open magnetic flux produced by a 
non-potential (NP) magnetic field evolution model for the 
lower solar corona (below 2.5iiQ). As observational input, 
this model uses the same observed photospheric magne- 
tograms as the widely used PFSS model, but while the PFSS 
open flux is determined entirely by the photospheric flux dis- 
tribution, the NP model predicts additional open magnetic 
flux due to inflation of the magnetic field by coronal electric 
currents. These currents arise naturally from the emergence 
of twisted active regions, their interaction, and their trans- 
port over the solar surface by photospheric motions. The 
total current varies over the solar cycle, leading to greater 
enhancement of the PFSS open flux at Maximum than at 
Minimum, as demonstrated in Figure 2. We propose that 
the presence of currents in the corona is a viable explana- 
tion for the activity-dependent discrepancy between PFSS 
predictions and in situ IMF observations at lAU. 

Of course, the influence of coronal currents as modeled 
here is only one possible explanation for the shortfall in 
PFSS open flux. To accurately determine the amount of 
open flux enhancement caused by coronal currents as op- 
posed, for example, to kinematic effects in the solar wind 
[Lockwood et al., 2009a], improved observational input is re- 
quired to better constrain the model. Firstly, discrepancies 
between PFSS extrapolations demonstrate already the need 
for accurate calibration of photospheric magnetographs (as 
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seen in Figure 1). Additionally, for the NP model, a key pa- 
rameter is the amount of electric current that emerges from 
the solar interior in bipolar active regions. In this paper, 
the twist of each newly-emerging bipole has been selected at 
random from a normal distribution approximating the ob- 
served overall distribution of active region helicity [Pevtsov 
et al., 1995]. Given this approximation, we have been able 
to compare only the global open flux with average IMF mea- 
surements. In particular, we have not considered the spatial 
distribution of open flux source regions {e.g., coronal holes) 
on the solar surface, which has potentially important impli- 
cations for acceleration of the solar wind [Cranmer, 2009]. 
To assess the more detailed impact of electric currents on 
the localised distribution of open flux will require reliable 
measurements of currents in individual active regions. Tech- 
niques both for extrapolating non-potential fields [DeRosa 
et al., 2009] and for estimating the rate of helicity injection 
[Demoulin and Pariat, 2009] from photospheric observations 
are under development. This effort will be greatly aided by 
the recently-launched NASA Solar Dynamics Observatory, 
which promises routine high-resolution, high-cadence, vec- 
tor magnetic field measurements at the photospheric level. 

A promising feature of the NP model is that it reduces the 
dependence of the open flux on the radius of the outer model 
boundary. In this paper, the boundary has been fixed at 
r = 2.5Rq, to allow comparison with PFSS extrapolations. 
However, this is an arbitrary choice made in the PFSS model 
to optimize agreement with observed coronal structures, and 
lacks physical motivation [Marsch and Richter, 1984; Jiang 
et al., 2010]. In the NP model, the magnetic field strength 
falls off less strongly with radius. We estimate on average 
that |Bo| ~ r~^'^ near the outer boundary, rather than 
as in the PFSS model, so that the open flux is nearer to be- 
ing independent of Rss- In future, we hope to remove the 
need for an artificial boundary altogether by coupling to a 
more realistic model of the solar wind. 

An important effect included in our NP model but not 
in existing extrapolation models is the ejection of twisted 
magnetic flux ropes, which form as a natural consequence 
of the transport of coronal magnetic field line footpoints by 
surface motions. Their ejection is a natural limiting mecha- 
nism for the build-up of currents in the corona [Bieber and 
Rust, 1995], acting on a timescale much faster than resis- 
tive dissipation. Moreover, it has been suggested that the 
temporary build-up of closed flux in the heliosphere from 
the resulting CMEs could be responsible for the discrep- 
ancy between PFSS extrapolations and IMF measurements 
[Riley, 2007; Owens et al, 2008]. Indeed, in our NP model, 
flux rope ejections cause fluctuations, on timescales of a few 
days, in the total open flux. However, our present quasi- 
static sinmlations with observational input only once per 27 
day Carrington rotation are unable to accurately estimate 
the timescale and level of enhancement to the Sun's open 
flux arising from individual CMEs. This will need to be 
considered in more detail in our future calibration of the 
NP model. For example, the additional contribution from 
ejections might counter-balance the reduction in open flux 
that would result from extending the outer boundary be- 
yond 2.5Rq. Ultimately, observations of CME contributions 
to the heliospheric magnetic flux could help to constrain the 
amount of electric current present in the low corona, and 
hence the most appropriate balance of relaxation and turbu- 
lent diffusion in the NP model. However, it is not clear that 
changing the timescales of these short-term fluctuations will 
strongly influence the background, mean level of total open 
flux in the NP model. This underlying enhancement of open 
flux arises from coronal currents absent in the PFSS model, 
and is a long-term phenomenon, reflecting the variation of 
the Sun's magnetic activity over months and years. 
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